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The S=l/2 Heisenberg chain with bond alternation and randomness of antiferromagnetic 
(AFM) and ferromagnetic (FM) interactions is investigated by quantum Monte Carlo simula- 
tions of loop/cluster algorithm. Our results have shown interesting finite temperature magnetic 
properties of this model. The relevance of our study to the observed results of the material 
(CH 3 )2CHNH 3 Cu(CLBri- ;c )3 is discussed. 
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I. INTRODUCTION 

Randomness induced quantum phase transitions have attracted intensive interests in the past decades. Putting 
enough strong bond randomness, analyses of real space renormalization group (RSRG) method have shown that RG 
flows of spin S=l and S=l/2 quantum antiferromagnetic Heisenberg chains 0,0, []| go to a stable fixed point called 
the random singlet (RS) phase in which spins far apart in space randomly form weakly bound singlet pairs. This 
property induces universal behaviors of ground states and low temperature thermodynamics, e.g. the energy spectrum 
is gapless, the temporal correlation length £ r and spatial correlation length diverge at zero temperature, and there 
is a non- universal infinite dynamical exponent z which comes from the relation ~ More important, the uniform 
susceptibility diverges universally in the RS phase as 

1 

y ~ fl) 

Tln 2 (ft/T) 

at low temperature, where f2 is a non-universal constant. On the other hand, the S=l/2 dimerized AFM chain was 
found to be extremely stable against strong bond randomness^. In this case, the system is in a quantum Griffiths- 
McCoy (QG) phase when the bond randomness is strong enough. This phase is characterized by gapless excitations 
and finite correlation length. In the QG phase, the uniform susceptibility behaves as 

X ~ T' 1 (2) 

at low temperature, where 7 is a non-universal exponent. On the experiments, the bond randomness effects 
have been found in several antiferromagnetic quasi-lD materials H H % Especially, a recent experiment on 
BaCu2(Sio. 5660.5)207 clearly show typical scaling relations of RS phase predicted by theory^. Moreover, there 
is another kind of bond randomness whose bonds can be both AFM and FM. In such systems, the RSRG analyses 
0, E| predicted a universal fixed point different from the RS phase because the spins correlate to form effective spins 
whose average size grows with lowering of the energy scale, and the magnetic susceptibility is Curie-like Xu ~ 1/T. 
Materials with such randomness is also fabricated 

Besides the above mentioned two kinds of bond randomness, Manaka and coworkers recently found that the com- 
pound (CH3)2CHNH3Cu(Cl x Bri_ :c )3 ^3] can be considered as a bond randomness S=l/2 AFM and FM alternating 
Heisenberg chains. The isomorphous compounds (CH3) 2 CHNH3CUCI3 13] is regarded as quasi-lD S=l/2 FM-AFM 
alternating material, and (CH3)2CHNH3CuBr3 14] is a S=l/2 AFM dimerized material. Mixing these two corn- 



ed 



pounds, they observed that a gapless phase appeared in the regime of the intermediate concentration 0.44 < x < 0.87 
of FM bonds by measuring magnetic susceptibility and specific heats. In order to describe the properties of this 
material, Hida 15] and Nakamura 0] suggested a ID model 
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where S, presents a spin S=l/2, J > 0, Jj = Jf(< 0) with a probability p — x 2 and Ji — Ja(> 0) with 1 — p. This 
model has two limits: (i) whenp = 1, it is a S=l/2 AFM-FM alternating Hciscnbcrg spin chain. When | J,/J |> 1 its 
ground state is the Haldane phase with gapped energy spectrum; (ii) when p — 0, it is a dimerized AFM Heisenberg 
spin chain, its ground state is the singlet dimer phase with gapped spectrum. By density matrix renormalization 
group (DMRG) method, Hida [13 considered the case p > 0.6 with J = 1.0 and | Ji |= 2.0 or 4.0, and confirmed that 
there exists QG singularity when p < 0.7. Nakamura [16J studied the model by non-equilibrium relaxation analysis of 
the quantum Monte Carlo (QMC) simulation in the whole parameter space of concentration < x < 1 with J = 1.0 
and | Ji |= 2.0, and found the gap vanishes in the regime 0.44 < x < 0.87 consistent with the experimental results 
|12| . However, in all these numerical works, the finite temperature magnetic properties in the whole parameter space 
of x are absent. In order to directly compared with the experimental results |12j . we perform the finite temperature 
QMC simulations extensively on this model in this paper. 



II. QMC SIMULATION RESULTS 

We investigate the magnetic and thermodynamic properties of the model defined by eq. with J = 1.0 and 
| Jj |= 2.0 using QMC simulations of continuous imaginary time loop/cluster algorithm. The weak AFM coupling 
J = 1.0 is fixed on all odd position bonds. For even position bonds, the strong AFM or FM bonds \ Ji \ = 2.0 are chosen 
randomly according to the probability p. We perform simulations for 100 ~ 200 bond arrangement configurations. 
For each bond configuration, after 2000 Monte Carlo steps (MCS) for thermalization, we further update 2000 MCS 
for Monte Carlo average. The temporal and spatial periodic boundary conditions are chosen for all simulations. In 
order to convince us of the code validity, we study S=l/2 dimerized AFM chain in weak dimerization, and find our 
results for the energy gap consist well with recent DMRG results 0] . 

We first investigate the ground state properties on the system size L = 128 for temperatures as low as (3 — 1/T = 200. 
At low temperature, the energy gap A is estimated by 



A = lim — , 

L — >oo F T 



(4) 



where £ r is the imaginary time correlation length obtained by second-moment method [T^- Then the valence-bond- 
solid (VBS) order parameter 0] 
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Z L =< exp[i— J^jS*] > 



(5) 



is measured to 



The results obtained are shown in Fig. ^ 
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FIG. 1: The energy gap A and FIG. 2: The illustration of spin or- 

VBS order parameter Zl versus dered phases with the probability 

the probability p. p — and 1.0. 

It is interesting to find that energy gap A exists for the Haldane phase at p — 1 point. It gradually approaches 
to zero at p w 0.7, and it opens again at very small p ks 0.02. It's hard to locate accurately the vanishing point of 
A because our A results are not exact zero value due to finite-size effects. But our calculations show that the gap 
should close at p c \ w 0.02 and p C 2 « 0.7. As a result, we presume that the system stays in the gapless phase in the 
regime of p c \ < p < p C 2 , which is consistent with the previous results 0] . 

In Fig. ^ we also observe that Zl ~ 1.0 at the limit p = 1, keeps at finite in the regime 0.35 < p < 1.0, changes its 
sign at p w 0.35, and then turns down to —1.0 when p < 0.35. The values of Zi « ±1.0 characterize the two limits of 
different ordered phases presented in Fig. [5] 
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Recent QMC study |20j on S=l random bond-alternating Haldane chain has shown that the VBS order parameter 
Zl is not effected by QG singularity, and it is an effective parameter to locate the RS critical point. Combining the 
results of Zl and A, we find that in the regime of 0.02 < p < 0.35 where A vanishes and Zl approaches zero from 
finite values, the system belongs to a critical phase. At the point p rs 0.35, as Zl changes its sign, a phase transition 
happens. In the regime of 0.35 < p < 0.7, A vanishes and Zi increases from zero to finite values. This fact reveals 
the system enters to other critical phase. In the regime of 0.7 < p < 1.0, where both A and Zl are finite, the system 
keeps in an ordered phase. 

In order to distinguish the upper mentioned different phases, we calculate the uniform magnetic susceptibility \u 
over the whole parameter space < p < 1.0 at finite temperatures. The results are summarized as following. 
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FIG. 3: The uniform susceptibility 
at 0.02 < p < 0.35 versus tempera- 
ture T = 1/13. 



FIG. 4: The fitness of Xu by 

ln~ 2 (fi/T) 



T — > 0, and every \u curve appears a valley, that is the typical feature of eq. for denoting RS phase. We fit the 
curves of p = 0.04 and p — 0.15 by eq. (JJJ, and find they can be very well fitted as plotted in Fig. 01 Thus we believe 
that the phase in this regime belongs to RS phase. 

If. The regime of 0.35 < p < 0.7 is also a gapless regime, where Xu curves diverge too, but they are obviously 
different from those in the regime of 0.02 < p < 0.35. Instead, as plotted in Fig. these curves are very similar 
as those in the ^v,o 0<3 +v, Q +-, r ™,-.ai f Q+ nl ^o io described bv ea. J2J. In FieEl we fit our low temperature 
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FIG. 5: The uniform susceptibility FIG. 6: The fitness of X u by T~ 7 

at 0.35 < p < 0.7 versus tempera- 
ture T = 1/(3. 

results of Xu by T -7 , and find again the fitness are quite good. It is interesting to note that the behavior of Xu is not 
Curie-like, so one can believe the phase is not belong to the universal class of AFM and FM bonds randomness H. 
We thus conjecture that the system is now in QG phase. 

III. For the regime of 0.7 < p < 1.0, the system enters to a gapped phase because all Xu curves appear the tendencies 
going to zero when T — > 0. Our results are plotted in Fig. [7| 

IV. At last, we consider the two limit cases p = and p = 1, and our results are plotted in Fig. [S] Obviously, the 
system belong to the gapped Haldane phase and dimerized AMF phase on these two points, respectively. 
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FIG. 7: The uniform susceptibility 
at 0.7 < p < 1.0 versus temperature 
T = 1/0. 




FIG. 8: The uniform susceptibility 
at p — and p = 1.0 versus tem- 
perature T = 1//3. 



III. CONCLUSION AND DISCUSSION 



From our QMC calculations, we can conclude that this S=l/2 AFM-FM alternating bond randomness chain has 
four different phases with respect to the probability p: (i) p = 0, the system is a dimcrizcd AFM chains with gapped 
energy spectrum; (ii) in the regime of 0.02 < p < 0.35, the system enters to the RS phase, whose energy spectrum is 
gapless and the uniform magnetic susceptibility Xu obeys the eq. {Q; (Hi) in the regime of 0.35 < p < 0.7, the system 
turns to the QG phase where the energy gap vanishes and the curves of Xu consist with eq. I|2"|): (iv) in the regime of 
0.7 < p < 1.0, the system is again in a gapped phase. Finally, the case p — 1 corresponds to the gapped AFM-FM 
alternating spin chain. 

Consequently, there should be three phase boundaries between these different phases: («)—>(«), because there is no 
effective quantities to locate the exact position of this boundary, we only say that the transition from the dimerized 
phase to the RS phase happens at very small p m 0.02; (ii)— *(iii), the phase boundary between RS phase and the QG 
phase resides at p = 0.35, where both the results of VBS order parameter Zl and susceptibility Xu consist reciprocally. 
Zl changing its sign at the point p = 0.35 also implies that it is a good quantity to indicate the transition from RS 
phase to other phase, which confirm the previous argument [20j; (Hi)— the location of the boundary between 
QG phase and the gapped phase at p w 0.7 is hard to be determined by results of energy gap A because of the 
finite-size-effects, but it can be extracted from the behaviors of Xu at different temperature. In the QG phase, Xu 
diverges as T — > 0, but in a gapped phase, Xu approach zero when T — > 0. Thus there should be a cross point of 
Xu at different temperature which correspond to the boundary. We investigate the case of same system size L = 128 
at different temperature 1/T — (3 = 10,50,100,200, and analyze the results by finite-size-scaling of imaginary time 
[2lj]. Our results, plotted in Fig. have i 1 .1.-- 1 ' • ' 0.71 ± 0.01, this result consist with 
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FIG. 9: Finite size scaling by T = 1/(3 of the uniform susceptibility near p c = 0.71. 
previous DMRG results [l5| . 

This model can be described under the scenario for the study of S=l AF bond random chain by R. Hyman and K. 
Yang [3J|, and our result give a positive evidence to support their conclusion. More important, our QMC results of Xu 
are consistent with the experiment results 0] of quasi-lD (CH3)2CHNH3Cu(Cl a; Bri_ a ;)3 in the two limits p = 0, 1, 
as found in Fig|U In the gapless regime, 0.02 < p < 0.71, where Xu diverge when T — > 0, the behavior are similar 
with the experiments observation that Xu(T) diverges also for 0.56 < x < 0.83 at low temperature. Especially, we 
find the critical point p c2 — 0.71 corresponding to x = 0.84 is very close to the experimental x = 0.83. As the critical 
point p c i w 0.02 corresponding to x ~ 0.14 is different from the experimental x = 0.56, it implies the current single 
chain model eq. © is failed to describe the results of experiment for small p, and some additional terms such as 
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the weak inter-chain coupling should be taken into account. Furthermore, we believe that the behavior of Xu{T) in 
the regime 0.56 < x < 0.87 for the experiment should successively exhibit first the QG-type divergence T~ 7 , then 
the RS-type divergence In -2 (Q/T)/T at low temperatures. Our further simulations are under consideration, and 
interesting results are expected in the near future. 
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